Bose-Einstein condensation in a decorated lattice: an application to supersolid 
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The Bose-Einstein condensation of vacancies in a three-dimensional decorated lattice is considered. 
The model describes possible scenario of superfluidity of solid helium, caused by the presence of 
zero-point vacancies in a dislocation network. It is shown that the temperature of Bose-Einstein 
condensation decreases under increase of the length of the segments of the network, and the law 
of decrease depends essentially on the properties of the vertexes of the network. If the vertexes 
correspond to barriers with a small transparency, the critical temperature is inversely as the square 
of the length of the segment. On the contrary, if the vertexes correspond to traps for the vacancies 
(it is energetically preferable for the vacancies to localize at the vertexes) , an exponential lowering of 
the temperature of transition takes place. The highest temperature of Bose-Einstein condensation is 
reached in the intermediate case of vertexes with large transparency, but in the absence of tendency 
of localization in them. In the latter case the critical temperature is inversely as the length of the 
segment. 

PACS numbers: 67.40.-w, 67.80.-s 

I. INTRODUCTION 

Experimental observation of non-classical rotational inertia in torsion experiments on solid He [l[ (confirmed by a 
number of other groups has revived interest to the idea on supersolid. The idea goes back to pioneer work by 

Andreev and Lifshitz [j| where it was shown that the presence of vacancies in quantum crystals at zero temperature 
(zero-point vacancies) can cause superfluid properties of such systems. However, as was found later [6, 7]. in He 
crystals the concentration of vacancies is negligible small and distincts from zero only due to thermal excitation. The 
presence of extensive defects in the crystal can change the situation and make the occurrence of zero-point vacancies 
energetically favorable. If an extensive defect is homogeneous, zero-point vacancies can move freely along the defect. 
At sufficient concentration of defects they form a network, that provides possibility of flowing of the vacancies through 
the whole crystal. Under lowering of the temperature such a gas of vacancies should go into a superfluid state. In 
recent papers [1, Q it was established by Monte Carlo simulation that in 4 He crystals the grain boundaries and 
dislocations do possess superfluid properties. On the other hand, as was shown in 0], a rather long-term annealing 
(that removes dislocation from the crystal) leads to a complete disappearance of the effect of a step-like change of the 
period of oscillation of a torsion pendulum filled with solid helium. Thus, superfluidity of vacancies in a network of 
dislocations can be considered as probable mechanism of superfluidit y o f quantum crystals. 

The idea on dislocation superfluidity was put forward in Refs. [lol. Illl] long before the observation of the effect 
As was shown in [ToL , the important parameter that determines the temperature of the phase transition in such a 
system is the length of the segment of the network. However the approach [l(| [H| did not consider such characteristic 
of the network as the transparency of the vertexes (more precisely, there was some implicit assumption on its value). 

In this paper we consider a simple model of the network that allows to investigate the dependence of the temperature 
of the Bose-Einstein condensation on the length of the segment and on tunnel characteristics of the vertexes. 

II. THE MODEL 

In what follows we will model a superfluid dislocation as one-dimensional lattice chain with the period a and the 
length I = qa. It is implied that the chain has nonzero concentration of zero-point vacancies n\ (one-dimensional 
concentration). The chains (segments) are joined into a regular three-dimensional network. Two edge sites of each 
segment are the vertexes of the network (we will call such sites the central ones). To be more specific, we consider 
that the network obeys cubic symmetry. In fact, we model the network as a decorated cubic lattice with the period 
I and the number of sites in the elementary cell equal 3(7 — 2 (Fig. [I}. We are interested in the case of large q. The 
vacancies that moves in a such lattice are described in the tight-binding approximation. The model contains two 
parameters: i, the amplitude of tunnelling between nearest neighbor internal sites in segments, and t\, the amplitude 
of tunnelling between a central site and a nearest neighbor internal site. The one-site energies are assumed to be the 
same for all sites. In the broader sense, the model describes a network formed by one-dimensional wires, along which 
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FIG. 1: Elementary cell of the decorated lattice (q = 5). Red circles represent central sites, green ones, the internal sites of the 
segments 



bosons can move. The central sites play the role of scatterers that connect the wires. 

The number of zero-point vacancies N is supposed to be much smaller than the numbers of sites in the decorated 
lattice, i.e. the filling factor for the vacancies satisfies the inequality v = n\a -C 1. For the further analysis it is 
convenient to introduce three-dimensional density of vacancies n ~ 3qis/l 3 . We note that even at small filling factor 
i/< 1 the number of vacancies per elementary cell can be much larger than unity (nl 3 3> 1) if the condition v ^> 1/q 
is satisfied. Below we will consider the filling factors belonging to the diapason l/5<f<l. 

At small filling factors one can neglect the interaction between the vacancies. The temperature of Bose-Einstein 
condensation Xq f° r a non-interacting gas of bosons in a compound lattice is determined by the equation 
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(1) 



where e.\(k) is the spectrum of bosons in a lattice, A, the band index, k, the wave vector, cq, the energy that 
corresponds to a bottom of the lowest band. In {1} the index A runs from 1 to 3g — 2 (the number of sites in the unit 
cell), and summation in taken over k belonging to the first Brillouin zone (— ir/l < fcj < w/l). 
The Hamiltonian of the system has the form 
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where bf n {bi yV ) is the operator of creation (annihilation) of a boson in the site r\ in the i-th cell, i, the radius- vector 
of the i-th cell, = (I, 0, 0), n y = (0, 1, 0), n 2 = (0, 0, 1), the primitive vectors of translation. The following notation 
for rj is used: rj = v, the central site, rj = (cc,£), the £-th internal internal site in a segment aligned in a direction. 
Applying the Fourier-transformation 
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(where iVj - number of unit cells), we rewrite the Hamiltonian as 
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The matrix M(k) has dimension (3g — 2) x (3g — 2) and is presented in the following block form: 
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Here D g _i is the (g — 1) x (g — 1) matrix that corresponds to the tunnelling between internal sites: 
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and T Q and T+ are of lx(q— 1) and (q — 1) x 1 matrixes that describe the tunnelling between the central site and 
the nearest internal site: 
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t = ii /t is a key parameter of the model (the ratio of amplitudes of tunnelling between an internal site and a central 
site and between two internal sites). 

The spectrum of bosons satisfies the dispersion equation det(el — M) = 0. Using the explicit expressions for M 
I©-®, we obtain the explicit form of the dispersion equation 

[A,_i(e)] 2 [eA,_i(e) - 6r 2 A 9 _ 2 (e~) - (-l)?2r 2 £ cos(fc Q l)] = 0, (8) 

a 

where e = e/t and A g (e) = det(el + D 9 ). 



III. BAND STRUCTURE OF THE SPECTRUM 



As follows from the dispersion equation (J8j> , some bands are reduced to degenerate levels. The energies of these 
levels are given by the equation 

A 9 _x(c) = 0. (9) 

Eq. |9]) coincides with the dispersion equation for an isolated chain with q — 1 sites. The energies of the levels are 
equal to 

7TS 

e s = -2tcos— (10) 

(s = 1, 2, . . . , q— 1). The degree of degeneracy of each level is 2iV,. The wave functions of such states has zero weight 
in the central sites that is in agreement with the absence of dispersion. 
The spectrum of the bands with a finite dispersion satisfies the equation 

eA g _r(e) - 6r 2 A 9 _ 2 (e~) - (-l)«2r 2 ^ cos(M) = 0. (11) 

a 

The left hand part of Eq. (jll[) is a q-th order polynomial. It has q distinct solutions that corresponds to q bands 
(allowed bands). The widths of the bands essentially depend on the parameter r. It is illustrated in Fig. [2l 

One can see from this figure, that there is a special value t — t c — 1/V3 at which all allowed bands have the 
largest width, and there is no gaps between the bands. At t > t c and t < t c the allowed bands becomes narrower 
and energy gaps (forbidden bands) appear. The effect of narrowing is strongest for the lowest and the highest bands. 
As we will see later, the temperature of Bose-Einstein condensation is determined in the main part by the width of 
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FIG. 2: The band structure for bosons in the decorated lattice (for q = 9). The allowed bands are shown grey. 



the lowest band. On a qualitative level the dependence of the width of bands on r is explained as follows. At small 
r it is energetically preferable for the vacancies to localize inside segments, and the central sites play the role of the 
barriers with a small transparency. As is known, in particular, from the Kronig-Penny model (see, for example, [l2l|). 
in such a situation a sharp narrowing of the lowest band takes place. At t > t c it is energetically preferable for the 
vacancies to localize at central sites, and overlapping between such localized states becomes very small. 

To obtain quantitative estimation for the Bosc-Einstein condensation temperatures it is necessary to find the 
spectrum of the bands for an arbitrary q. From the definition of A g (e) follows that this function satisfies the recurrent 
relation 

A q (e) = eA q . 1 (e)-A q . 2 (e) (12) 

(A^e) = e, A 2 (e) = e 2 — 1). Using the relation (TT2"]) and applying the method of a mathematical induction one can 
prove that 

A g (2cos 7 ) = Sin[(g+1)7] (13) 
sin 7 

and 
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Eqs. (fT"3f and (fT4|) allow to reduce the equations (|9|) and (fTT|) to compact trigonometrical equations. In particular, 
the substitution 1=2 cos 7 reduces Eq. ^ to the equation sin 57 = (with the additional condition sin 7^0) which 
solutions correspond to the energies (fT0|) . 

For finding the spectrum we use the substitution e = 2 cos(7r + 7) that reduces Eq. (|ll|l to the following equation 
for 7: 



cot 7 sin(g7)(l - 3r 2 ) + 3r 2 ^cos(<77) - jjj cos(k a l)j 



0. (15) 



At t = t c = l/y/3 Eq. USD can be easily solves that yields the spectrum 
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(j = 1, 2, ... q), where square brackets indicate the integer part. One can see from (|16l) that at any q the bandgaps 
are equal to zero. We display also approximate expression for the half-width of lower bands 
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(j C?). As follows from (fT7|) . the half- widths are in inverse proportion with the square of length of the segment (at 
fixed a). 

At r < t c one can find an approximate solutions of Eq. (|15[) in a diapason of small 7 (that correspond to the lower 
bands). We specify the case of r not too close to r c , when the condition (1 — 3r 2 )g S> 1 is satisfied. The solution of 
Eq. (|15p can sought as 7 = 7^ + jj, where jj = irj j q (j -C q) and 7j <C jj. The approximate expression for jj as 
obtained from (fT"5l) reads as 



3r 2 



7j 



1 - 



(-I)' 



1 - 3r 2 q 2 v 

It gives the following expression for the spectrum of lower bands 
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In this case the half-width of the lowest band is equal to 
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i.e., it is inversely proportional the third power of the length of the segment. We note that in the Kronig-Penny 
model at small transparency of the barriers the width of the bottom band is also inversely proportional to the cube 
of distance between the barriers (see, for instance, [13]). 

If r exceeds r c , (and g>l)a sharp narrowing of the lowest band takes place, and this band drops below the level 
—2t. Using the substitution e = — 2 cosh 7, we obtain from (fTTj) the following equation 



coth7 sinh(#7)(l - 3r 2 ) + 3t 2 
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(21) 



At q(3r 2 — 1) ^> 1 it is convenient to seek for a solution of Eq. (|2Tj) in the form 7 = 70 + 7, where 70 is given by the 
equation 



coth7o = 



3r 2 



3r 2 - 1 



(i.e. 70 = ln(6r 2 — l)/2), and an exponentially small correction 7 can be obtained directly from (|2"Tj) 

7 ~ - sinh7o COSI170 e~ qla cos(k a l). 



(22) 



(23) 



The spectrum of the lowest band reads as 



ei(k)w-£b 1 



(24) 



where 



E = 2t cosh 70 = 2t 



3r 2 



V6r 2 - 1 



Thus, the lowest band becomes exponentially narrow (with the half-width W\ « At cosh 70 sinh 2 7 e^« 7 °). All the 
others bands lay above the level — 2i. The spectra for the lower bands (starting from j = 2) are determined by the 
equation Q19p if one replaces j for j + 1 in its left hand part. The eigenfunctions for the lowest band have the maximum 
weight in the central site and it falls quickly with the distance from the central site. 
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IV. TEMPERATURE OF BOSE- EINSTEIN CONDENSATION 



Due to the presence of degenerate levels in the spectrum the formula |T]) for the temperature of Bose-Einstein 
condensation(BEC) is modified to 

9 1 1 

N = YY t— , +2NiY 7 n ■ (25) 

fe^exp(^)_l feex P (W)-l 

We are interested in systems with a rather large concentration of vacancies (n ^ l~ 3 ). In this case the BEC 
temperature is much larger than the width of the lowest band. It allows to use the approximate expression for the 
Bose distribution function for such a band 

1 T 



exp( 6l( ^~ £0 ) -1 ei( k )- £ o' 



The rest bands and degenerate levels give non- negligible contribution into ([25]) , only if their energies, counted from 
a bottom of the lowest band, are of order or less than the BEC temperature. The main contribution into (j2"5)) yields 
the lowest band. Therefore, for evaluation of what bands and levels should be taken into account, one can use the 
estimate T ~ Winl 3 . 

Let us first consider r < r c (and q ^> 1/(1 — t 2 /t 2 )). In this range of the parameters the half-width of the lowest 
band W\ cx q~ 3 , while the energy gap between the first and second bands A g ^ cx q~ 2 . At small filling factors v the 
three-dimensional density satisfies the condition nl 3 = Zqv <C q. Hence, To <C A g .i and it is enough to take into 
account in Eq. ([25]) only the lowest band and the lowest degenerate level (that lies at the top of the lowest band) . As 
a result, Eq. (f2"5)) is reduced to 

n/ 3 «— / / / dk x dk v dk z : = — = = — + — «2.5— . (26) 

Wo Jo Jo iyi[l-i(cosfc x +cosfc y + cosfc,)] Wx Wi 

It is instructive to write down the answer for the BEC temperature in terms of one-dimensional concentration of 
vacancies m, the distance between nearest sites a, and the length of the segment I: 

To = 1.2^^2*^. (27) 

At r rs t c the energy gaps between the bands approach zero and one should take into account many bands and 
levels. Since the main contribution into (|25p is given by the levels and the band, which energy counted from the 
bottom of the lowest band is smaller than To, it is enough to take into account only such levels and bands. Besides, 
it is possible to approximate all bands, except the first one, by the degenerate levels (with degree of degeneracy A/j), 
located between degenerate levels (|TU| . As a result, Eq. (|2"5T) is reduced to 

„/3_ J o 1 / / / ,7 ,7 ,7 I . v- la . in 



nl -WF^ / / dk x dk v dk z — +2 > — u —+^ ^—3, (28) 

W^ 3 Jo Jo Jo 2(larccos^7 fc °) ^ 2W ' S U^{j + \f 

where s max and j m ax - numbers of levels and bands, which energy is of order of To. In view of fast convergence of the 
sums in (|28[) they can be extended to infinity. The integral in Eq. (f2"B")) can be evaluated numerically. As a result, we 
obtain T « 0.2Winl 3 . Using Q2J> we find 

T w3ta 2 y (29) 

At r > t c the energy gap that separates the lowest band from the lowest level and the second band is exponentially 
large in comparison with the width of the bottom band. Hence, it is enough to consider in (I25p only the lowest band. 
It gives T w 0.7Winl 3 . At large q = l/a this temperature is exponentially small. For example, at r = 1 the BEC 
temperature is equal to 

T w 9We~^ (30) 

We should note that since the states in the lowest band correspond to the vacancies localized near the central site, 
the result (f3"0")) is more sensitive to the interaction than results for other r. Due to the tendency to localization at 
central sites, the interaction can be neglected, only if number of vacancies per elementary cell is less or of order of 
unity. Therefore, the result (|3T))) should be considered as qualitative one. 
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V. CONCLUSION 



The model of decorated lattice considered in this paper describes three physically distinct situations, depending on 
value of t. The case of small r corresponds the situation, when intersections of dislocations play the role of barriers 
with a small transparency for the vacancies. In this case the BEC temperature is inversely proportional to the square 
length of the segment. The case r > r c describes the situation, where the vacancies tend to localize at the dislocation 
intersections. In such a situation the BEC temperature decreases exponentially with the increase of the length of the 
segment. At last, the case r *=s r c corresponds the situation where vacancies can move freely through the intersections, 
and do not localize on them. Last situation is the most preferable for the BEC: only linear decrease of the temperature 
of transition with the increase of the length of the segment takes place. We consider that the functional dependences 
of the BEC temperature on the length of the segment is more or less universal and does not depend on the mechanism 
that causes the appearance of the barriers (or the centers of localization) on the dislocation intersections. 

It is of interest to estimate BEC temperature for the most preferable case (121)1) . If one defines the effective mass of 
vacancies on a simple (not decorated) cubic lattice M* = h 2 /2ta 2 the result (|29j) can be presented in the form 

3 h 2 m 

T ° w 2iiFT (31) 

Note that the answer (pJTj) does not contain any parameter of the decorated lattice. Thus this answer can be applied 
to a more general case, where dislocations can be modelled as quasi one-dimensional wires in which vacancies can 
move. The answer ([31]) up to the numerical factor of order of unity coincides with the estimate given in Ref. [111 ]. 
Implying that the effective mass of vacancies is approximately equal to the mass of 4 He atom, we evaluate (|3"TT) as 
To w (ni/l) ■ 18 K -A 2 . For example, at density of vacancies n\ = 1 A -1 (a value obtained in 0]) the temperature 
To w 0.1 K is reached for the length of the segment / = 180 A, that corresponds to two-dimensional density of 
dislocations rid = 3 • 10 11 cm~ 2 . 

In conclusion, we discuss shortly the question on the relation between the Bose-Einstein condensation and the 
superfluidity of vacancies in the network of dislocations. At temperatures lower than the BEC temperature the 
long-range phase correlations are established in the system. In our case the phases on different segments become 
correlated. It allows to describe the system (below To) by a complex order parameter that reflects the possibility of 
non-dissipative flow along the dislocations. Above To the flowing without relaxation is impossible. But the specifics of 
the system studied consists in that the temperature To and the temperature of degeneracy for the gas of vacancies on 
a segment are quite different: Td = tn\a 2 >• T . Therefore, the relaxation time for the flow at ^> T > T ; can 
be very large, and above T a peculiar quasi-superfluid phase can be realized (see [III) and also discussion in 0). Let 
us note that in the experiment (l3j the attempt of direct observation of Bose-Einstein condensation in solid helium 
was made (through the measurements of temperature dependence of the kinetic energy per atom). According to the 
results of Ref. [13J, in the temperature diapason where nonclassical rotational inertia is observed (measurements were 
done down to 0.07 K), the Bose-Einstein condensate does not arise yet. 
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